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ABSTRACT 

We have analysed a large data set of VI absorber candidates found in the 
spectra of 3702 SDSS quasars, focusing on a subsample of 387 AGN sightlines 
with an average S/N> 5.0, allowing for the detection of absorbers above a rest- 
frame equivalent width limit of W r > 0.19A for the O VI 1032 A component. 
Accounting for random interlopers mimicking an O VI doublet, we derive for the 
first time a secure lower limit for the redshift number density AN/Az for red- 
shifts Zaiyg > 2.8. With extensive Monte Carlo simulations we quantify the losses 
of absorbers due to blending with the ubiquitous Lya forest lines, and estimate 
the success rate of retrieving each individual candidate as a function of its red- 
shift, the emission redshift of the quasar, the strength of the absorber and the 
measured S/N of the spectrum by modelling typical Ly forest spectra. These cor- 
rection factors allow us to derive the 'incompleteness and S/N corrected' redshift 
number densities of O VI absorbers :ANovi,c/ 'Az c (2.8 < z < 3.2) = 4.6 ± 0.3, 
AN OV i,c/Az c {3.2 < z < 3.6) = 6.7 ± 0.8,and AA ov/jC /Az c (3.6 < z < 4.0) = 
8.4 ±2.9. 

We can place a secure lower limit for the contribution of O VI to the closure 
mass density at the redshifts probed here: Qovi(2-8 < z < 3.2) > 1.9 x 10~ 8 /i _1 . 
We show that the strong lines we probe account for over 65% of the mass in the 
O VI absorbers; the weak absorbers, while dominant in line number density, do 
not contribute significantly to the mass density. 
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Making a conservative assumption about the ionisation fraction, and adopt- 



er 



ing the lAnders fc Grevessd (119891 ) solar abundance values, we derive the mean 
metallicty of the gas probed in our search : ((2.8 < z < 3.2) > 3.6 x 10~ 4 h, in 
good agreement with other studies. These results demonstrate that large spec- 
troscopic datasets such as SDSS can play an important role in QSO absorption 
line studies, in spite of the relatively low resolution. 



Subject headings: Quasar Absorption Lines - IGM 



Introduction 



The cosmic evolution of the majority of the baryons, as probed by absorption features 
of neutral hydrogen towards the sightlines of QSOs, has come to be well understood over 
the past years of extensive research : the availability of powerful echelle spectrographs on 
8m-class telescopes like HIRES/Keck and UVES/VLT has opened the wind ow for high res- 



olution studies of the Lyman a forest down to very low column densities ( iHu et al.l Il995 



Lu et al.lll996l ; iKim et al.l 119971 ; iKirkman fc Tytlerl 119971 ). On the theoretical front, models 
incorporating gas dynamics, radiative cooling, and photoionisation have been spectacularly 
successful in reproducing the observed properties of the parts of the gas traced by neutral 



HI throughout much of the history of the universe flCen et a 



1994; Miralda-Escude et al. 



19961 ; iHernquist et al.lll996l ; iPetitjean et al.lll995l ; iDave et al.lll997f ). The emerging picture 
paints baryonic gas being encountered in a wide variety of physical conditions. Most of 
the high redshift Lyman a forest, which is the main repository for baryons at that time of 
the universe, exhibits relatively cool temperatures (T~ 10 4 K) in low-to-medium overden- 
sity structures, governed by photoionisation, but there are clear signs of a additional highly 
ionised component, traced by the prominent absorption lines of CIV and OVI. 
The existence of metal absorption lines in the spectra of quasars , already well established 
early on in the study of QSO absorption lines ( iBahcall et al.lll968[ ). and then later linked to 
the fo rest (C IV absorbers seen at the same redshift as Lyman a absorbers by lMeyer fc York 
( 119871 ). poses interesting problems, and in fact allows to address several astrophysical phe- 
nomena in a unique setting. First of all, metals allow insight into the structure of the 
intergalactic medium (IGM) itself : its temperature and small-scale velocity structures can 
be probed via the width of such lines. Secondly, the existence of metal absorbers reveals 
the impact of star formation on the IGM - either the metals had to be produced insitu 
by Population III stars in small-scale regions, or a transport mechanism from galaxies like 
galactic winds into the IGM has to be invoked. Furthermore, studying the ratios of metal 
column densities of different species and ionic stages can be used to constrain the spectral 
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shape of the metagalactic UV background. 

The hot phase of the IGM (WHIM) with temperatures above T ~ 10 5 K is still compara- 
tively poorly understood. Because collisional ionisation of hydrogen becomes significant at 
these temperatures, Lyman a becomes a worse tracer of this phase, and thus budgets of the 
IGM content based solely on hydrogen will underestimate the contribution of the hot gas 
(ISimcoe et al.ll2004j ). The high ionisation potential and high abundance renders oxygen a 
prime candidate for studying the IGM phase under these conditions, possibly tracing shock- 
heated material. Furthermore, OVI has been long recognised as ideal candidate to trace the 
low density component (log N#j < 14.0) of the photo ionised Lyman forest at z> 2 with 
optical ground-based spectroscopy (IHellsten et al.lll998l ). with C IV being another possible 
tracer. The low density in this regime allows for the p roduction of OVI via ph otoionisation 
by the intergalactic UV and soft X-ray radiation field (IHaardt fc Madaulll996l ). 
Searches for the signatures of intervening OVI absorbers in QSO spectra have focused on 
two techniques. A direct line detection, based upon the identification of the characteristic 
OVI doublet at rest-frame wavelengths of 1032 and 1038 A, at low redshifts must be per- 
formed via space-based observations, and ground-based detections hence become possible 
only for absorber redshifts b eyond z n h.« ~ 2.0 , when the lines are shifted redward of the 
atmospheric cutoff at 3000 A. iBurles fc Tytlerl f 19961 ) conducted the first systematic survey 
for OVI absorbers at z a f, s ~1.0 with FOS aboard HST, studying a sample of 11 QSOs and 
identifying 9 good candidates, thus establishing that the number density per unit redshift 
of these absorbers is similar to, i f not greater than, tha t for CIV and Mgll absorbers at 
the same redshift. More recently, iThom fc Chen! (j2008al Jbl) have used HST/STIS to search 
for intervening OVI absorbers towards 16 QSOs at z em = 0.17 — 0.57, retrieving a sample 
of 27 candidates with rest-frame equivalent widths W r >5mA. The b iggest sample of such 
absorbers at low z a b s has been collected by iDanforth fc Shulll (120071 ). who obtain 83 OVI 
systems by analysing 650 Lya absorbers with HST and FUSE spectroscopy, enabling sound 
statistical analyses of the population. 

Towards higher redshift, despite the advantage of having the features shifted into the wave- 
bands for ground-based optical spectroscopy, this search method is hampered by the ever 
increasing density of the Lya forest, rendering an unambiguous detection of the absorber 
doublet in its midst more difficult. Hence, most searches have employed high-resolution and 
high signal-to-noise (S/N) spectroscopy of a handful of QSO sightlines, to build up samples 
of such OVI absorbers. In order to ensure the highest likelihood of detecting the doublet 
securely, these studies have focused on absorber and emitter redshifts below z ~3.0. To date, 



range 2.0 < z n h. < 3.0 is of order 350 ( 


Simcoe et al. 


2004 




2006; 


Bergeron & Herbert-Fort 


2005; 


Carswell et al. 


2002; 


Fox et al. 


2007; 


Lopez et al 


. 2007 


Fox et al. 


2008 





An alternative approach to characterising the metal content of the Lyman forest is the 
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pixel-to-pixel-correlation method, pioneered by ICowie fc Songailal ( 119981 ) for CIV. With this 
method, statistical correlations between the mean optical depths of various transitions en- 
able an estimate of the metallicity of low density g as. Over the s a me redshift range as 
mentioned for the direct line identifi cation searches, ISongailal (119981 ); ISchaye et al.l (120001 ) 
and most recently lAracil et al.l ( 120051 ) assert the existence of OVI in gas of optical depth as 
low as thi = 0.1. For a summary of th e res ults and a discuss i on of the ambiguity in their 
interpretations see lAguirre et al.l (120021) andlPieri fc Haehneltl (120041 ) . 

In the first paper of our series (IFrank et al.ll2010l ). we have described how we utilise the 
large number of high-redshift QSOs in the Sloan Digital Sky Survey (SDSS) spectroscopic 
sample to explore a different locus of the parameter space for direct searches : low resolution 
(R ~ 1800) and low to medium S/N at emission redshifts z em > 2.8. Frank et al. (2010) 
describes the sample selection, the search algorithm and presents the candidates found in 
over 3700 sightlines. While it is clear that the lower quality of the spectra as compared to 
the high-resolution studies does not allow us to probe such weak absorbers as in the above 
mentioned searches, the large combined pathlength covered by the SDSS QSOs' sightlines 
enables us to systematically build up a large sample of strong O VI absorbers at unprece- 
dented high redshifts. 

In this paper, we explore the statistics of the subsample of the candidates towards sightlines 
with the best S/N, and derive, for the first time at redshift beyond z a b s > 2.8, secure lower 
limits for the redshift path density, dn/dz, of strong O VI absorbers. Furthermore, by a 
Monte-Carlo analysis we can determine the fraction of absorbers lost due to noise and the 
blending with ubiquitious Lyman forest lines, and hence derive a good estimate for the path 
density probed with our analysis (section 2). Enforcing a strict lower limit to the total col- 
umn density found for our sample, we derive in section 3 a limit for the total mass density of 
O VI at these high redshifts, and then in section 4 we estimate the minimum metallicity of 
the gas traced by our absorber sample by limiting the ionisation of oxygen to the maximum 
value achieved in both the collisionally and photoionised cases. We compare our results with 
other studies, and conclude in section 5, where we also propose different scenarios for the 
origin of these absorbers and point out the need for higher resolution and/or signal-to- noise 
spectra in order to discriminate between them. 

Unless specifically noted otherwise, we use a cosmology with H =71 km s -1 Mpc -1 , Qm — 
0.27 and Q\= 0.73 in this study. Abunda nces are given by number relative to hydrogen, and 
we utilise the lAnders fc Grevessd ( 119891 ) solar abundance values. 
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2. Frequency of O VI absorbers 



2.1. The uncorrected Redshift Number Density 

The redshift number density of VI absorbers dn/dz, defined as the number of absorbers 
per unit redshift, is a function of the limiting equivalent width W r j imit . Here we estimate 
dn/dz as : 

n to t{> W r) u mit ) - n rand (> W r 

limit ) i , \ 



dnldz{W r ,umit) 



Az t , 



ot 



where n tnt {> H ',..;,,„,/) i s the total number of absorbers found by the search procedure de- 

0|), n randi> ^r,iimit) is the estimated number of random interlop- 
and Aztot the total pathlength covered by all spectra used for 



scribed in Frank et al. 



ers and false classifications 
assembling the sample. The error on this redshift density is 



O 'dn/dz 



V^rand 



Az t , 



(2) 



ot 



The total pathlength for a sample can be found by summing up all sightlines that cover 
various intervals of redshift. In our case, all AGN spectra that were searched for VI 
absorbers had the same lo w redshift boundary, defined by the instrumental setup of the SDSS 
spectroscopy ( IYorkll2000l ). In order to avoid including intrinsic and associated absorbers into 
our sample, we have chosen to define the maximum redshift such that the velocity difference 
between the AGN's emission and the absorber is 5000 km/s (where we define blue-shifted 
absorbers to have positive velocity values). This 'arbitrary' velocity cut-off proves to select 
against most of the galaxies clustered around the QSO (and ionized by the QSO) ; all the gas 



in th e > 1 kpc region of the QSO host, and most of the "intrinsic" absorbers (jWild et al. 
20081 ). But note that it has become clear that a small fraction of the "intrinsic" absorbers 



can have (3 up to 12,000 km/s (IVanden Berk et all ( 120081 ) and iNestor et all fj2008h . Hence, 
limiting ourselves to the lower cut-off of 5000 km/s will, in the strict sense, provide only 
upper limits for the number of 'intervening' absorbers above this cut-oflE, but we estimate 
that the errors introduced by the unknown fraction of 'intrinsic' absorbers with (5 > 5000 
km/s are negligible compared to the other systematic and statistical errors of our algorithm. 
Fig. 1 shows g(z), the number of lines-of-sight within z and z + dz, from which the redshift 
path Az can be easily obtained : 

Az = Y>g(z)5z (3) 



1 In order to estimate the number of interlopers and the fraction of true lines lost in the forest, we have 
created a large set of artificial spectra. For details see next section. 

2 Note that the opposite case, intervening absorbers at (3 <5000 km/s is also possible. Hence, the true 
total number density of intervening absorbers cannot be assessed here. 
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where 5z is the resolution of the grid. Note that the overwhelming majority of sight-lines in 
the full sample is comprised of spectra with a signal-to-noise ratio (SNR) within the Lyman 
forest below the nominal value of 4.0, that an AGN needs to obtain overall to enter the SDSS 
QSO database. In order to ensure comparability with other studies, and to allow for the 
unambiguous detection of O VI absorption at a reasonable significance level, we have chosen 
to use a subsample of the full dataset which includes only sources for which SNR/ orest > 5.0. 
In the following, this is going to be the sample from which we draw inferences unless otherwise 
stated. It consists of 387 AGN, and allows for the detection of O VI absorbers above a rest- 
frame equivalent width limit of the 1032 A component of W r > 0.19A. 

Table 1 lists the number of AGN, the total pathlength, the total number of absorbers 
found, as well as the estimated number of random absorbers per redshift interval, and the 
line density derived from these quantities. For comparison it also lists the findings of various 
other studies for O VI and other absorbers. Note that these values represent the observed line 
densities based upon actual detection of lines, no correction for the fraction assumed to have 
gone undetected has been applied. Thus, in a strict sense the measurements only provide 
a secure lower limit. For the full sample over the complete redshift range 2.8 < z a t, s < 5.1, 
i.e. including all sources regardless of S/N within the forest , a density An/Az = 0.21 ± 0.01 



is obtained. This is much lower than the value derived by iBurles &: Tytlerl (119961 ) for their 
sample at an average redshift of z a h s = 0.9, but this can be explained both by the fact 
that our identifica tion efficiency degrades rapidly towards higher redshifts (as detailed in 
Frank et al.l (120101 ). and that, as noted above, there does not exist a clearly quantifiable 
equivalent width limit W r> u m it for the majority of the sources where SNR^ orest < 4.0. When 
focusing on possible absorbers in the redshift interval 2.8 < z a {, s < 3.2 within the sample 
of high SNR/ ores t, where we expect the detection efficiency to be best, our line density of 
A N/ Az(2.8 < z„h{,3.2) = 1.31 ± 0.12 is in good agreement with the lower redshift value 
of IBurles fc Tytlerl (119961 ). indicating that there is no evolution of the rate of incidence of 
str ong O VI absorbers ove r this wide range of redshifts. We stress, however, that although 
the lBurles &: Tytlerl (119961 ) analysis was performed with similar spectral resolution and SNR, 
the fact that the Lyman a forest density is much higher at the redshifts probed in our work, 
probably leads to a much more severe underestimate of the line density in our case due to 
the higher fraction of lines lost in the forest by blending with strong H I absorbers, as we 
will see in the next section. 

The comparison to searches for O VI absorbers at only slightly lower redshifts than probed 
by our study is difficult, since most of these were performed using high SNR and much 
higher spe ctral resolution than the SDSS spectra allowed . Thus, it is not surprising that 
studies like Isimcoe et all (120041 ) ; iBergeron fc Herbert-Fortl d2005h ; ICarswell et aD (booj ) find 
a substantially higher An/Az at z a b s ~ 2.0 — 2.5: the limiting equivalent width in their 
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Fig. I. — The number of sightlines towards AGN in the full SDSS sample, and the subset of 
sources where the SNRj orest > 5.0, versus the maximum redshift allowed in the spectrum for 
detecting intervening absorbers. We have chosen a velocity difference between the absorber's 
redshift and the quasar emission of (5c > 5000 km/s. 
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cases drops to W ri u m it ~ 0.05A, allowing to probe the more abundant weaker systems. 
Simply dropping their weak absorbers from their samples to derive the frequency of strong 
O VI absorbers does not provide a reliable estimate, because multiple weak and for SDSS 
unresolvable components of a complex absorbing system can form a feature detectable with 
our search algorithm, yet the number of lines identified in such cases in the high resoultion 
studies strongly depends on the often detailed velocity and column density structures of the 
components. Hence, in a strict sense, we need to compare the number of 'systems', and not 
the individual components of these studies with our 'can didates' that all represent strong 
features. Of more importance, however, is the fact that iBergeron Sz Herbert-Fortl (120051 ) 
employ a different method in estimating An/Az, they infer the line density by integrating 
over a fit to the column density distribution f(N): An/Az = AX/ Az x J f(N)dN, where 
AX is the redshift path density. This procedure requires detailed knowledge of the column 
density sampled by each absorber, and can usually only be obtained with high resolution 
and high S/N, as we will detail in the next section. 



2.2. The incompleteness and S/N corrected Redshift Number Density 



As indicated in the preceding section, the 'raw' number density obtained by simply 
counting the O VI candidates found with our search algorithm has to be modified in order 
to find the 'true' number of such absorbers. As a first step, we have already incorporated 
an estimate of the number of spurious H I interlopers that mimick O VI doublets. This 
estimate stems from running our search algorithm on all spectra with a wavelength separation 
between the two putative absorption features that has a random offset from t he no minal 
rest-frame AA = 1037.613A - 1031. 912A = 5.70lA, as detailed in iFrank et all (bold ). Far 
more important, however, is the estimate of fraction of candidates lost through blends with 
ubiquitious Lyman forest lines. For each candidate we have calculated this fraction by 
creating a set of mock forest spectra, as we describe below. 

We take the absorber redshift, the equivalent width estimate, the AGN emission redshift and 
the S/N me asurement inside t he relevant part of the Lyman forest from our catalogue list of 
candidates ({Frank et al.l 120101 ). Then we create 200 differe nt line lists , rand omly populated 
with H I absorbers according t o the line d ensity estimate of iKim et al.l (120011 ) with a column 
density distribution following iHu et al.l (119951 ) . for the appropriate emission redshift z em 
covering the complete wavelength range from the SDSS low-cutoff (~ 3800 A) up to the 
Lyman a emission of the quasar. From this list we create a mock high-resolution spectrum 
by modelling the Lyman absorbers with Voigt profiles, where the Doppler parameter b is 
drawn from a Gaussian distribution with a mean of b c = 25.0 km/s and width of A = 8 
km/s. We have modelled all H I transitions down to Ly e. We then degrade the spectrum 
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Table 1. Observed Line Density Of Various Absorption Line Studies 



Ion 


Zabs 


Sight-lines 


r, limit 


Aztot 


n toi 




An/Az 


Refer( 


VI 


2.8 - 3.2 


305 


0.19 A 


98.9 


139 


9.8 


1.31 ± 0.12 


1 


VI 


3.2 - 3.6 


84 


0.19 A 


53.2 


41 


3.5 


0.71 ± 0.14 


1 


VI 


3.6 - 4.0 


9 


0.19 A 


9.74 


5 


0.63 


0.45 ± 0.24 


1 


VI 


2.8- 5.1 


3366 




3768.4 


1290 


508.9 


0.21 ± 0.01 


2 


VI 


0.9 


11 


0.21 A 


4.73 


6 


1.3 


1.0 ± 0.6 


3 


VI 


2.2 - 2.5 


7 


0.04 A 


2.016 


230 * 


n.a. 


114 


4 


VI 


2.2 - 2.5 


7 


0.04 A 


2.016 


-100 ** 


n.a. 


50 ± 5 


4 


VI 


1.99 - 2.57 


10 


? A 




136 


n.a. 


74 (66.. 106) 


5 


VI 


1.99 - 2.57 


10 


0.05 A 


? 


? 


n.a. 


26 


5 


VI 


< 0.5 


16 


0.03 A 


2.623 


55 


n.a. 


21.0±3.0 


6 


VI 


< 0.5 


16 


0.30 A 


3.179 


1 


n.a. 


0.3 +0.7/-0.3 


6 


VI 


2.1-2.3 


2 


0.03 A 


0.70 


8 


n.a. 


11.4 ± 4.0 


7 



References. - Note that lines density, as used by the authors, refers more strictly to 
components of systems, which are all counted as separate entities. 

1. This work, SNR>5 sample, 

2. This work, full sample , W ri u m it not constant., 
3. 
4. 



Buries fc Tvtlerl (119961) . 



Simcoe et al.l (120041 ) - high SNR + Resolution, slightly diffferent method (survival analy- 



sis ), * including also upper l i mits, ** including detected lines only 

5. iBergeron fc Herbert-Fortl J2005f ) - high SNR + Resolution, from fit to CDDF, 

no EW limit given for first value, in second row no values given for number of lines and 
pa thlength 

6. iTripp et al.l (120081 ) , no correction for false identifications 7. ICarswell et al.l ( 120021 ). 
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to the SDSS resolution, add the appropriate level of noise to each pixel, and place the O VI 
candidate into this mode l of the Lyman for est. Running the same detection algorithm on 



these spectra as used in iFrank et al.l (120101 ). we can thus estimate how often the Lyman 
forest destroys the signal. The fraction fj of absorbers retrievable is therefore a function of 
the absorber redshift, the emission redshift of the quasar, the average S/N within the forest, 
and the strength of the absorber as measured by its equivalent width : 

fi = fiizabs, z em} S/N } W r (1032A), T^(1038A)) (4) 

We can then proceed to estimate the 'true' number n to t,c of absorbers that we should have 
found by summing over all candidates i : 

n to t,c( W r > W rMmit ) = j (5) 

Exactly the same procedure can be performed for the estimate of the number of random 
interlopers : 

n ra nd,c( W r > W r ,Umit) = T ^ 

o h 

where the summation is over all j candidates found via the 'random' search. 
Additionally, we have taken into account that there are individual pixels in the spectra that 
will fall below the S/N limit we have imposed in the first step for the average S/N inside 
the Lyman forest. We thus have to reduce each individual redshift path Azi towards the 
ith quasar by removing such 'bad' pixels. This results in a lower overall Az to t,c for each 
redshift interval of interest. This effect is much more pronounced towards higher redshifts, 
as the rapidly increasing density of the Lyman forest introduces more and more regions with 
substantial flux decrements, leading to lower individual S/N at the affected wavelengths. 
It is interesting to see that these corrections reverse the trend within the redshift bins avail- 
able to us of a decreasing number density with increasing redshift : while we now arrive at 
An c /Az c {2.8 < z abs < 3.2) = 4.6 ± 0.3, the values for 3.2 < z abs < 3.6, and 3.6 < z abs < 4.0 
increase to An c /Az c = 6.7 ±0.8, and 8.4 ±2.9, respectively. As Table 2 shows this is mainly 
due to the vastly decreased corrected redshift path length. Due to the large errors, and the 
risk of being susceptible to cosmic variance with such short pathlengths, we cannot assess 
the reality of this trend. If, however, this trend was real, then it would imply an extremely 
strong evolution of the redshift path density : An c /Az c ~ (1 ± ;z) 7 , where 7 ~ 2.5 — 3.0, 
compar able or even strong e r than th e evolution of e. g . stro ng Mg II absorbers at z abs < 2.3 



seen by ISteidel fc Sargentl (119921 ) or iProchter et al.l (120061 ). A comparison to the studies 



mentioned in the preceding section is not possible, because none of these has incorporated 
all the corrections we have applied here, and also probably sample different absorbers, as we 
indicate in the next section. 
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Fig. 2. — The OVI line density estimates for our and other studies. Our data are repre- 
sented by the three redshifts bins above z a ^ s > 2.8. Indicated are the secure lower limits 
on An/Az that we can place simply by directly counting our candidates and summing up 
all redshift paths (arrows). By accounting for the losses of absorbers in the dense forest via 
analysing a large set of mock spectra, and the reduced pathlength induced by the noise in 
certain pixels being higher than allowed for our detection limit of W r > 0.19A, we are able to 
derive the corrected values of An/Az shown for the same bins (points with error bars). The 
measurements and/or upper and lower limits for other studies are labeled, and their redshift 
range is indicated by the length of the bins. Note that the three studies at intermediate red- 
shifts (2. < z < 2.7) are perfor med with high-resolution spectroscopy (HIRES or UVES), 
while the iBurles fc Tytler! (119961 ) is done with similar resolution as ours (FOS/HST). The 
upper limit of dn/dz = 114 indicated for the Simcoe sample (blue dotted lines) 
includes all of their upper limits/non-detections. The data point for Bergeron 
(red dashed) is their conservative estimate when applying a rest-frame equiva- 
lent width limit of 50mA, dropping this limit they obtain dn/dz=74(+32-6), as 
indicated by the dashed upper limit. For details see text. 
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3. The Cosmological Mass Density 

The contribution of a certain trace element in a particular ionisation state to the cur- 
rent closure density of the universe can be estimated from the measurement of the column 
densities (Lanzetta et al. 1991, Songaila et al. 2001). In our case : 

O - 1 m S N OVIj ( 7 \ 

where p c = 1.89 x 10~ 29 g cm -3 is the current critical density, movi the mass of the O VI 
ion, Hq the current Hubble parameter, and Novi,i the total O VI column density toward the 
zth QSO. The absorption path distance AX between the maximum and minimum redshifts 
of a sight-line towards a QSO, [z min , z max ], for a flat universe is defined as 

AX = [Q m (l + z max f + n x } 1 ^ 2 - [tt m (l + z mm f + fi A ] 1/2 (8) 

The standard error of this mass density can be estimated to : 



= 1 y/ZiWovi,- < Nqv7>? 

au [1 - (1/n)] - 5 AX tot { ) 

where n is the total number of absorption systems in the sample, and < Nqvi > is the mean 
column density of the sample (Buries 1996). 

A lower limit on the column density Novi i n each absorbing system can be obtained 
by measuring the rest-frame equivalent width of the weaker 1038 A transition W r (1038 A). 
Fig. [3] shows the curve-of-growth for this ion under various model assumptions. If the 
absorption features are unsaturated, the column density and rest-frame equivalent width are 
simply proportional to each other. For saturated lines, however, the relation retains a model 
dependency : larger broadening parameters b lead to larger W r . In a first step, we have 
extracted iy r (1038 A) from the catalogue of all lines identified by our pixel-by-pixel search 
algorithm (Frank et al. (2010) in the high SNR sample of spectra (SNR/ orest >5.0, 387 
spectra, 237 O VI identifications). In order to decide whether two absorption features at the 

Table 2. Incompleteness and S/N corrected Line Density Estimates 



Ion 


Zabs 


Sight-lines 


r, limit 




n tot,c 


H-rand,c 


An c / Az c 


O VI 


2.8 - 3.2 


305 


0.19 A 


50.5 


251.7 


19.4 


4.6 ± 0.3 


O VI 


3.2 - 3.6 


84 


0.19 A 


12.3 


92.7 


10.3 


6.7 ± 0.8 


O VI 


3.6 - 4.0 


9 


0.19 A 


1.37 


13.75 


2.2 


8.4 ± 2.9 
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Fig. 3. — The curve of growth for OVI 1038 A. The relation between the rest-frame equivalent 
width W r and the column density Nqvi is obtained by modelling the 1038 A transition with a 
Voigt profile with a given Doppler parameter b. Note that the linear part the curve of growth 
does not depend on the choice of b, whereas the flat part retains that model dependency. 
The solid lines on the left side show the 25th percentile, average, and 75th 
percentile, while the dotted lines indicate the extreme values of the equivalent 
width distribution of the O VI lines in the SNR>5 sample. 
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correct wavelength separation can qualify as VI 1032/1038 A doublet, the automatic search 
program 'measures' the equivalent width of an absorption feature simply by multiplying the 
flux decrement in one pixel with the pixel wavelength separation. This gives, as detailed in 
Frank et al. (2007), a lower limit to the equivalent width as an absorption feature might 
span more than one pixel. In Fig. [3] the location of the average W r (1038 A) and the extreme 
values found in our sample are denoted by the solid and dotted lines. It is obvious that even 
these lower limits for the equivalent widths begin to populate the regime where the reliable 
estimate of the column density from W r alone is not possible. We can, however, draw sound 
infer ences on lower limit s : Fig. [3] shows that even for unreasonably high b paramaters (note 
that ISimcoe et all (120021 ) obtain a median b = 16 km/s) very few of the absorbers can reside 
below Novi < 10 14 cm -2 , and the average for the sample is well above that limit. Thus, we 
can safely assume that ^ Novi,i > i x 10 14 cm" 2 . 

Table [3] shows the results of these calculations. For the redshift bin with the highest 
confidence and the best success rate of the automatic search algorithm, we obtain a lower 
limit for the mass density of fiov/(2.8 < z a bs < 3.2) > 1.9 x 10 _8 /i _1 . How does this 
value compare to measurements at lower redshifts ? Buries & Tytler (1996), analysing FOS 
spectra with similar resolution, arrive at an O VI density limit at an average redshift of 
z ave = 0.9 of Qovi > 9. 8 ±3.0 x 10 8 /i 1 (but note that they use a different cosmology). The 
inspection of two QSO sightlines with high resolution UVES spectra by Carswell, Schaye 
and Kim (2002) yields Qovi = 6 x 10~ 8 /i _1 at z ~ 2.2, and the analysis of Simcoe et 
al. ( 2004) for a sample of z ah .« ~ 2.5 abs orbers gives Vt OVI = 7.89 x 10~ 8 x h" 1 . Note 
that iBergeron &: Herbert-Fortl ((2005) give a much higher mass density over the 
range 2.0 < z a b s < 2.6, when accounting for incompleteness via integrating over 
the column density distribution function : VLovi = 24.5 x 10~ 8 x h~ l . Hence, there 
is a wide range of values found by different authors using different methods at 
slightly lower redshifts than probed by us. Apart from the Buries & Tytler (1996) 
study, these values are derived from analyses of high-resolution and high-S/N data. And yet 
our estimate based upon SDSS spectra already allows us to pose a lower limit at a higher 
redshift that lies only a factor of 3-5 below the mass density esimates of Carswell, 
Schaye and Kim (2002), and the one of Simcoe et al. (2004). It is interesting 
to note that the incompleteness correction terms introduced for our line density 
estimate are of that orde r. We remain, however, a factor of 13 lower than 



Bergeron fc Herbert-Fortl ( 120051 ). and hence whether there is a trend suggesting 



for the O VI mass density to decrease when moving from z a b s ~ 2.5 to z a b s ~ 3.2, is 
thus not completely clear. This also means that the more abundant, yet weaker 
O VI absorbers probed by the Simcoe et al. (2004) study cannot contribute 
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much to the O VI mass density^ 

Indeed, a closer examination of the cumulative fractions of the line density and the cumulative 
contribution to the baryon density of the absorbers found in our SDSS search reveals that 
we sample already a high fraction of the mass, despite missing a large fraction of the lines as 
we cannot detect weak absorbers. Figure [3] shows the fraction of VI absorbers detected as 
a function of the lower limit to the column density sampled (upper panel). This cumulative 
fraction, and the contribution to the baryon density of these absorbers (lower panel), 
are derived assuming a power-la w dependence for the column density distribution function 
( [Bergeron &: Herbert-Fortl 120051 ) : fN OVI ~ ^ovr The nne density per redshift interval, 
dn/dz, is then obtained by simply integrating the column density function : 



dn/dz{N low < Nov i < N, 



up J 



f(N)dN ~ N' 



-a+l 



(10) 



where the upper boundary for the O VI column density, N up , for the absorbers can be 
estimated from theoretical models for the IGM or assumed from published lists, and the 
lower boundary, Ni ow , is in principle a function of the detectability limit, varying with 
instrumental set-up and the redshift intervals probed. H ere we have chosen to adopt the 
slope a = 1.71, found by lBergeron &: Herbert-Fortl (120051 ) for redshifts z a b s ~ 2.5, assuming 
that there is no substantial redshift evolution. For the upper limit to the column density, 
we apply two different values. We have checked the published lists of O VI absorbers, and 
found that there are only very few examples of lines as strong as log Novi > 15.5. Even for 
damped-Lyman a (DLA) absorbers with their high column density for neutra l hydr ogen, the 
O VI absorbers associated with them tend to fall below this limit iFox et al.l ( 120061 ). Hence, 
we limit ourselves in the first case to column de nsities below Ngyj = 1 155 . A second, 



sensible limit can be derived by integrating the iBergeron fc Herbert-Fortl (120051 ) column 
density distribution, and calculating the number of strong absorbers we expect to be able 
to retrieve within our given pathlength. Beyond column densities of log Novi — 16.0, this 
number drops below three absorbers in the redshift range 2.8 < z a t, s < 3.2, and below one 
at higher redshift. We have thus chosen log Nqvi — 16.0 as our limiting column density 
for the second case. We stress that these upper limits probably reflect a physical boundary 
: when the gas density increases beyond a peak value for each species, the recombination 
drives the ionisation equilibrium fraction for the O VI ionic species to lower values than 
can be compensated for by the increasing gas density, and hence the O VI density will also 



3 Note that these authors also specifically make an incompleteness correction by integrating over a column 
density distribution function and a metallicity function, that was obtained by a survival analysis taking into 
account non-detections. 
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drop (in other words : recombination is dependent on n^, whereas ionisation only depends 
linearly on n e )Q There is no physical limit for the lower boundary, but lines too weak simply 
become undetectable. A column density of logN OVI = 12.0 produces an absorption feature 
that is right at the detection limit of typical high- resolution and high S/N exposures with 
8m-class telescopes and instruments like UVES or HIRES. Hence, we have set this value as 
our lower boundary Ni ow . Note that the lower boundary is much more important for the 
fraction of lines that can be discovered, than for the cumulative mass fraction sampled, as 
the power- law index a = 1.71 is close to 2. 

The baryonic mass traced by all detected lines within the interval given by the two boundaries 
discussed above, is proportional to the intergral 

N up 

tt b (N low < Novi < N up ) ~ J Nf(N)dN~N- a+2 (11) 

-^low 

Figure |3] shows in the upper panel the cumulative fraction of the line density as a function 
of the limiting column density in a survey. The red dotted line and arrow indicate how low 
our SDSS based sample reaches, and the inter val indicated by the blue dotted lines and the 



arrows delineates the boundaries chosen by ( iBergeron fc Herbert-Fortl 120051 ) to maximise 



completeness and minimise the effects of cosmic variance in their sample of UVES spectra. 
Note that our long effective pathlength allows us in principle to probe out to higher column 
density systems without being afflicted by the rarity of such absorbers causing problems 
with cosmic variance. Clearly, we cannot expect to sample more than 5% of the lines that 
can be retrieved with typical high- resolution and high-S/N observations. As the lower panel 
shows, however, we are picking up the bulk of the mass density of O VI absorbers with our 
approach : depending on the exact upper limit to the column density of such absorbers, we 
are sensitive to at least 65% of the absorbers by mass, and assuming the same lower cut-off 
as for the number density (log Nov/ > 12.0), we are in principle able to sample up to 
82% of the mass, not accounting for the inevitable incompleteness of the survey 
towards the lower end of the column density thresholdjf] 



4 Note that this is not true for the case of collisional ionisation. 

5 While the fraction of lines that can be retrieved is a strong function of the integration limits for the column 
density, the fraction of the mass sampled depends mostly on the high-cutoff value. But even restricting 
ourselves to the regime between 13.0 < logNovi < 15.0 means that we are probing a high fraction of the 
mass, as the magenta curves in this figure demonstrate. 
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Fig. 4. — The cumulative fractions of the line density dn/dz (upper panel) and the con- 



distribution for the column density f(N) ~ N 



Bereeron & Herbert-Fort ( 


2005), a = 1.71. 


ts of our and the 


Bergeron & Herbert-Fort 



( 120051 ) analyses, respectively. Note that while we do not retrieve a substantial fraction of all 
lines, we nevertheless sample > 65% of the mass in such absorbers. For details see text. 
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4. A lower limit for the Cosmic Metallicity 



From the cosmological mass density at redshift z, it is a straightforward proces s to 
estimate the mean metallicity of the gas due to OVI. Following iBurles fc Tytlerl (ll996f ). let 
C(z) be the ratio of mean metallicitjo at redshift z to the solar value, then 



^ {Z> ~ 16 { )Solar OVI 



;i2) 



where /x is the mean molecular weight (~ 1.3 for a 75% hydrogen and 25% helium mix by 



mass), and (^)soiar = 1174.9 (jAnders fc Grevessd Il989l ). The cosmological baryon mass 



density as measured by (ITytler fc Fan! 1 19941 ) from high-z deuterium is fib = 0.023/i -2 . With 
the current data, the ionisation fraction -^j cannot be determined, but we are able to place 
a str ict upper limit for th e O VI fraction of 22%, both for collisionally excited or photoionised 
gas (jAguirre et al.l 120071 ). 



With the lower limit for the O VI mass density from above, we obtain 

C(2.8 <z< 3.2) > 3.6 x 10~ 4 /i 
or expressed in the usual logarithmic units, 



(13) 



[O/H] > -3.6 



(14) 



Buries and Tytler (1996 and erratum 2002) derive ((z abs ~ 0.9) > 4.0 ± 1.2 x 10" 4 x h , 
but note that these authors operate with an ionization fraction of OVI/O = 1. 
Adjusting to our our correction factor, their result (C(z a bs ~ 0.9) > 18.2 ± 5.5 x 
10 -4 x h) is a factor of 5 higher than our (not very constraining) lower limit.. 
Carswell, Schaye and Kim (2002) obtain ((z a b. s ~ 2.2) ~ 16 x 10~ 4 x h, whereas Simcoe 
et al. (2004) propose a median [O/H] = -2.82 (i.e. ((z ~ 2.5) = 22 x 10~ 4 x h) for the 
cosmological filaments enriched with oxygen and carbon, but also see examples of systems 
with very low abundances. They estimate that up to 30% of the Lyman forest lines have 
[O/H] < -3.5, i.e. ( < 4.5 x 10" 4 x h. Using th e novel Loc a lly C alibrated Pixel 
(LCP) search method in >2000 SDSS QSOs, IPieri et al.l (EoS) find that the 
results of the search for OVI allows for tightly constraining the mass-weighted 
mean oxygen abundance for the IGM ([<0/H> w ] = -2.33 ± 0.23 (2a), while 
the volume-weighted mean abundance is less well constrained : [<0/H>yvy] = -2.74 



6 Note that this definition takes into account only the gas containing strong OVI systems accompanied by 
HI. In addition, it implicitly assumes a volume filling factor of one for the gas with such detected systems in 
SDSS spectra. Given these caveats, it is still a safe lower limit that we derive. 
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± 0.64 (2a). Note, however, that this technique is specifically tailored to the low 
density regime of the IGM, whereas our approach is complementary in the sense 
that we are dealing with strong absorbers arising presumably in much denser 
structures than Pieri et al. (2010) probe. 



5. Summary and Conclusions 



We have analysed the iFrank et al.l (120101 ) data set of O VI absorber candidates found 



in the spectra of 3800 SDSS quasars, focusing on a subsample of 387 AGN with an average 
S/N> 5.0 per pixel inside the Lyman forest region, allowing for the detection of absorbers 
above a rest-frame equivalent width limit of W r > 0.19A for the O VI 1032 A component. 
Accounting for random interlopers mimicking an O VI doublet, we derive for the first time 
for redshifts z a bs > 2.8 a secure lower limit for the redshift number density An/Az : 
An OV i/&z(2.% < z < 3.2) > 1.31 
An OVI /Az(i.2 < z < 3.6) > 0.71 
An OV i/Az{3.6 < z < 4.0) > 0.45 

In order to quantify the losses of absorbers due to blending with the ubiquitous Lya forest 
lines, we have modelled typical Lyman forests following canonical line density and column 
distribution functions, and estimate the success rate of retrieving each individual candidate 
as a function of its redshift, the emission redshift of the quasar, the strength of the absorber 
and the measured S/N for the detection spectrum. Furthermore, we have excluded the parts 
of the redshift pathlength in each spectrum that fall below the S/N which allows for the 
detection of lines above the rest-frame W r = 0.19A. These correction factors allow us to 
derive the 'incompleteness and S/N corrected' redshift number densities of O VI absorbers : 
An v/, c /Az c (2.8 < z < 3.2) = 4.6 ± 0.3 
An v/, c /Az c (3.2 < z < 3.6) = 6.7 ± 0.8 
An OV i,c/Az c {3.6 < z < 4.0) = 8.4 ± 2.9 

These values are lower than the ones found in other studies at redshifts 2 < z a b s < 2.5, owing 
to the fact that the low SDSS resolution does not allow us to probe systems as weak as in 
those high- resolution and high S/N studies. Simply comparing these studies with our ap- 
proach by dropping the weak lines out of the statistics for the high-resolution is not possible 
: multiple weak and thus for SDSS unresolvable components of a complex absorber can form 
features that will be detected by our algorithm, yet the number of lines identified in such 
cases strongly depends on the details of the velocity and column density structures. 
Transforming the measured equivalent width of the O VI candidates into a column density 
is not possible with the limited information obtainable with the SDSS resolution. We can, 
however, place a secure lower limit for the detectable O VI column density : log Nqvi > 14.0. 
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With this limit, we can put a constraint on the contribution of O VI to the closure mass 
density at the redshifts probed here: 
n OV i(2.8 < z < 3.2) > 1.9 x lO- 8 /*- 1 

Interestingly, this limit is only a factor of 2-3 below the density estimates of the high- 
resolution studies mentioned above. Noting that the correction factors due to blending with 
Lja and due to the decrease of the available path because of noise are roughly of the same 
order, this indicates that the weak O VI absorbers do not contribute significantly to the 
mass density. Indeed, we have shown that the strong O VI systems that we probe account 
for > 65% of fib(0 VI). The large number of sightlines allows us for detection of strong O VI 
absorbers despite the low number- density of the lines, while also overcoming the problem 
of cosmic variance. While high-resolution, high S/N studies of QSO absorption lines are 
invaluable for IGM studies, we have shown here that low resolution studies also provide 
important and complementary results, as long as they cover a large number of sightlines. 
Making the conservative assumption that the ionisation fraction, can never be above 
0.22 (for both the photoionised and the collisionally ionised cases), and adopting the Anders 
& Grevesse solar abundance values, we derive for the mean metallicity of the gas probed in 
our search : 

C(2.8 < z < 3.2) > 3.6 x 10" 4 /i 

in good agreement with other studies. We take this fact, i.e. that the limits on number den- 
sity, mass density and, to a lesser degree, the mean gas metallicity we derive agree 
well with these studies, as further indication that our search strategy yields sensible 
results. 

Having obtained these observational data will allow us to confront them with simulations 
that model the intergalactic medium's metal enrichment, and serve as important constraint 
for some of the model parameters. We point out that a similar strategy can easily be applied 
to similar, even larger data-sets (as e.g. the BOSS survey project for SDSSIII will provide 
in the near future), which will enable without doubt much smaller (statistical) errors. 

The SDSS is managed by the Astrophysical Research Consortium for the Participating 
Institutions. The Participating Institutions are the American Museum of Natural History, 
Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western 
Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Ad- 
vanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute 
for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the 
Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National 
Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for 
Astrophysics (MPA), New Mexico State University, Ohio State University, University of 
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Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Obser- 
vatory, and the University of Washington. 
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Table 3. Cosmological Mass Density of Various Absorption Line Studies 



Ion 


Zabs 


Sight-lines 


AX tot 


x h- 1 x 10- 


8 


Reference 


VI 


2.8 - 3.2 


305 


113.9 


> 1.9 




1 


VI 


3.2 - 3.6 


84 


24.1 


> 2.6 




1 


VI 


3.6 - 4.0 


9 


2.6 


> 2.9 




1 


VI 


0.9 


11 


6.6 


> 14.0± 4.0 (q = 


0.5) 


2 


VI 


0.9 


11 


6.6 


> 9.8± 3.0 (q = 


0.0) 


2 


VI 


2.5 


7 


? 


7.89 




3 


C IV 


2.5 


7 


3.27 


1.78 




3 


VI 


2.3 


10 


? 


15.1 




4 


VI 


2.3 


10 


? 


35(26..67) 




4 (fit) 



References. - 1. This work, SNR>5 sample, 2. Buries & Tytler 1996 - note different 
cosmology !, 3. Simcoe 2004, using Survival method. 4. Bergeron & Herbert-Fort (2005), 
using a fit to the column density distribution, and two different sample criteria 



